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Monitoring  Declining  Quality  of  Ammunition  Stockpile  under 

Step-Stress 

So  Young  Sohn 
Dept,  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Abstract 

Most  ammunition  is  produced  long  before  its  ultimate  consumption  and 
stored  in  a  series  of  different  depots  for  a  considerably  long  period  of  time. 
During  storage,  the  quality  of  the  ammunition  stockpile  deteriorates  pro¬ 
portionally  to  the  conditions  of  depots.  We  view  different  conditions  as¬ 
sociated  with  a  series  of  depots  as  step-stress.  A  random  effects  logistic 
regression  model  is  employed  to  predict  the  quality  of  ammunition  stock¬ 
pile  in  terms  of  the  routing  information  such  as  a  series  of  location  and 
duration  of  storage  of  ammunition  lots.  The  resultant  prediction  model 
can  be  used  to  determine  the  appropriate  time  for  reorder  or  reno\’ation  of 
ammunition  before  the  quality  reaches  substandard.  An  example  is  given  to 
illustrate  the  implementation  procedure  of  the  prediction  model  suggested 
in  this  paper. 

Key  Words:  Random  Effects  Model,  Step-Stress,  Deterioration  Rate, 
A  Two-Stage  Estimation 
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1.  INTRODUCTION 

'  Designing  the  proper  surveillance  program  for  the  material  whose  qual¬ 

ity  deteriorates  during  storage  has  been  one  of  the  important  topics  in  the 
area  of  material  management  (Valdez-Flores  and  Feldman  [13],  Whitehead 
[14]).  Especially,  when  the  degradation  of  the  quality  would  cause  not  on¬ 
ly  economic  loss  but  also  catastrophic  disaster  such  as  loss  of  human  life, 
importance  of  the  appropriate  quality  control  cannot  be  overemphasized. 
A  good  example  would  be  ammunition  lots  that  are  stored  in  depots  for  a 
relatively  long  period  of  time  before  their  ultimate  usage.  In  order  to  keep 
ammunition  stockpile  from  reaching  substandard  faster  than  expected,  well 
planned  surveillance  programs  are  necessary. 

In  an  attempt  to  provide  inputs  to  such  surveillance  programs  Sohn 
[10]  formulated  an  ammunition  stockpile  deterioration  model  based  on  a 
random  effect  logistic  regression  analysis.  The  suggested  model  enables 

*  one  to  predict  the  deterioration  rate  of  ammunition  lots  in  terms  of  depot 
condition  along  with  other  related  characteristics  such  as  vendor  sources 
and  the  manufacturing  year.  In  order  to  model  deteriorating  patterns,  one 
of  the  assumptions  employed  in  Sohn  [10]  was  that  once  ammunition  lots 
are  sent  to  a  depot,  they  remain  in  the  same  depot  during  the  experimental 
period. 

In  practice,  however,  ammunition  lots  may  be  transferred  to  several 
depots  in  sequence  before  ultimate  usage.  A  typical  logistics  chain  of  am¬ 
munition  lots  described  in  Brzuskiewicz  and  Morrison  [1]  begins  with  the 
load  plant.  At  the  load  plant,  ammunition  lots  sometimes  spend  up  to 
one  year  in  temporary  storage  before  deployment  to  the  permanent  storage 

•  area.  The  ammunition  lots  shipped  to  and  stored  in  permanent  depots  are 
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then  rotated  from  war  reserve  and  tested  on  schedules  which  depend  on  the 
availability  of  the  item  and  the  policy  of  the  Defense  Ammunition  Director. 

In  sum,  an  ammunition  lot  would  be  exposed  to  a  series  of  different 
conditions  of  depots  associated  with  different  levels  of  average  temperature 
and  humidity  in  various  locations  (underground  depot,  aboveground  depot, 
and  warship)(see  Eriksen  and  Str0msoe[2],  Forsyth  et  al.  [3]).  A  series  of 
different  conditions  can  be  viewed  as  step-stress  consisting  of  several  levels 
of  constant  stress. 

In  this  paper,  we  formulate  the  deteriorating  pattern  of  the  quality 
of  ammunition  stockpile  under  step-stress  in  order  to  accommodate  the 
possible  exposure  of  ammunition  lots  to  several  different  depots.  The  main 
goal  is  to  provide  tools  to  predict  the  quality  of  ammunition  given  the 
expected  duration  and  locations  of  storage  of  ammunition  lots  in  sequence. 
In  section  2,  the  model  formulation  (a  random  effects  logistic  regression) 
under  constant  stress  and  the  necessary  estimation  methods  introduced  in 
(Sohn  [10])  is  briefly  reviewed.  In  section  3,  by  combining  the  segments  of 
individual  models,  a  prediction  model  for  ammunition  deterioration  under 
step-stress  is  derived.  In  section  4,  an  example  is  given  to  illustrate  the 
implementation  procedure.  Finally,  discussion  is  given  in  section  5. 

2.  CONSTANT  STRESS  MODEL 

Consider  the  following  experiment.  N  lots  of  homogeneous  caliber  am¬ 
munition  (say,  fuze  manufactured  by  the  same  vendor  in  the  same  year)  are 
piurchased  from  a  vendor.  They  are  sent  to  m  different  depots;  lot  number 
1, ..,  ATj  to  depot  1;  lot  number  N^  +  l, N1+N2  to  depot  2;  ,..,  ;  lot  number 
N„,-^  1  +  1,  A^n.- 1  +  Nfn  =  N  to  depot  m. 

It  is  assumed  that  once  lots  are  stored  in  a  depot,  they  remain  in  the 
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same  depot  diiring  the  experimental  period.  As  a  result  of  acceptance 
sampling,  the  qualities  of  incoming  ammunition  lots  are  assumed  to  be 
homogeneous  while  the  average  deterioration  rates  of  ammunition  lots  may 
vary  depending  upon  conditions  of  depot  in  which  ammunition  lots  are 
stored.  Although  environmental  conditions  associated  with  a  depot  vary 
continually  over  time,  the  average  condition  of  a  depot  is  considered  as 
constant  stress  given  to  ammunition  lots  in  the  depot. 

In  order  to  estimate  deteriorating  patterns  of  ammunition  lots  over  time 
under  constant  stress,  each  lot  i  (i  =  l,.,iV)  is  repeatedly  inspected  on  a 
sampling  basis  without  rectification.  The  number  of  defective  items  ) 
found  out  of  sample  size  (riy)  at  the  jth  inspection  of  lot  i  (j  =  1, ..,  7i,) 
would  follow  a  binomial  distribution  with  a  parameter,  expected  cumulative 
proportion  defective  (py).  The  expected  cumulative  proportion  defective 
(Pij)  would  be  a  non-decreasing  function  of  time  (ty)  and  we  use  the  fol¬ 
lowing  within-lot  logistic  model  to  describe  the  deteriorating  pattern  of  lot 
i: 

For  f  =  1, ..,  fV,  and  j  —  I, ..,  7i, 

Pij  =  exp(J3o  +  0ity)/{l  -I-  expifla  -I- 0it^))  (1) 

where  exp(/3o)/(l  +  exp{l3o))  represents  the  initial  proportion  defective  of 
ammunition  lots  and  /3,  is  the  deterioration  rate  of  ammunition  lot  i  which 
would  be  positive.  The  average  deterioration  rate  of  ammunition  lots  stored 
in  one  depot  often  differs  from  that  of  another  depot  depending  upon  their 
environn'.ental  conditions.  We  assume  it  is  mainly  due  to  different  depot 
conditions  while  there  is  some  part  of  variation  that  cannot  be  explained 
by  such  conditions. 
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One  of  the  possible  models  that  accommodate  these  points  would  be 
the  following  between-lot  model  for  /?,: 

/3,  =  exp{')iZii  +  ..  +  7m2im  +  Ci)  or 


In^i  =  7i2.i  +  ..  +  7Vn2.„.  +  e,  (2) 

where  za,  {k  =  1, m)  is  a  dummy  variable  (2*=!  if  a  lot  i  is  stored 
in  depot  k]  otherwise  0)  and  e,  follows  independent  N(0,  cr^).  The  corre¬ 
sponding  regression  coefficient  7^  would  imply  the  average  log(deterioration 
rates)  of  ammunition  lots  stored  in  depot  k. 

Once  unknown  regression  fficients  (7'i,..,7m)  in  (2)  are  estimated, 
they  can  be  used  as  a  basis  to  examine  which  depots  are  associated  with 
significantly  higher  average  deterioration  rates  than  the  others  as  shown 
in  [9].  In  addition,  when  f3o  is  available,  one  can  predict  the  proportion 
defective  of  an  ammunition  lot  which  would  be  stored  in  one  of  the  depots 
or  similar  depots  used  in  the  experiment. 

In  order  to  estimate  unknown  parameters  such  as  /?o  and  (71, 7^),  we 
use  a  two-stage  method  which  separates  the  estimation  of  the  within-lot 
model  (1)  from  that  of  the  between-lot  model  (2)  (Korn  and  Whittemore 
[5],  Sohn  [10]  and  Stiratelli  et  al.  [11]). 

A  Two-Stage  Estimation 

First,  in  order  to  estimate  the  within-individual  model  (1),  the  following 
likelihood  function  of  i/,j  conditional  on  f3o,  is  formulated: 

ii=nii  (3) 

t=l  j~l 

where  is  as  in  equation  (1). 
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In  order  to  estimate  /3o,  I3i,...,I3n  that  maximizes  (3),  we  differentiate 
ln{Ll)  with  respect  to  (3q,  By  solving  a  set  of  resulting  normal 

equations,  we  find  the  maximum  likelihood  (ML)  estimates,  (/3o,  /3i,...,/3a-). 
A  (iV  + 1)  X  1  vector  4—  (/ioi  0i,---,0Ny  would  follow  asymptotically  normal 
distribution  with  mean  /3  =  (/3o,  and  variance  A  where 

0  0 

''T.O  "fi,!  ^1,2 


A  = 

^  '^Nfi  .  '’n,N  > 

The  inverse  of  the  negative  information  matrix  evaluated  at  (/?o, 

A  can  be  used  to  estimate  the  variance  matrix  A. 

Once  the  ML  estimates  4i’s  are  obtained,  they  can  replace  the  unobserv¬ 
able  Pi  in  the  between-individual  model  (2).  This  replacement,  however, 
adds  the  estimation  error  6,  to  the  equation  (2): 


In^J,  "f"  •••“!”  ^m^im  “b  (4) 

for  i  =  1, TV. 

Using  matrix  notation  model  (4)  can  be  written  as  follows: 

In/?  =  Z'y+  e-t-  (5  (5) 

where  ln0  is  an  iV  x  1  vector,  {hipi,  ..,iTV?.v)';  Z  is  an  N  x  m  matrix  of 
ZikS]  the  7  is  an  m  X  1  vector  of  7’s;  e  is  an  TV  x  1  vector  of  e,’s;  and  6  is 
an  iV  X  1  vector  of  (5,’s. 
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6  is  assumed  to  be  statistically  independent  of  e  and  would  asymptoti¬ 
cally  follow  normal  distribution  with  mean  0  and  variance  Q.  where 


Q  = 


^'^a/PnI3i 


'T2,nI02Pn 


^S,N 


/0N^  / 


Estimated  i?,  Q,  can  be  obtained  by  replacing  with  ^^k/Pi3k 

for  i,  k  =  1, N.  In  sum,  In/?  ~  N{Z'y,  V)  where  V=  Q  +  a-. 

Based  on  normal  In ,3,  the  restricted  maximum  likelihood  (REML)  method 
is  used  to  estimate  the  between-individual  model  parameter  (t~  which  sat¬ 
isfies  the  following: 


-  0.5  trR  -  O.SlnB'RR  InB  =  0  (6) 

where  R  =  -  V-^Z(Z' V-^Z)-  ^Z' V-^  (Searle  et.  al  [9]).  Once  d-  is 

obtained  from  (6),  it  replaces  cr^  in  V and  7  can  be  estimated  as 

7  =  {Z'  WZ)-\Z'  lV(/n/?))  (7) 

where  W=  {Q-\-  and  the  estimated  variance  of  7  is  (Z'  WZ)~^. 

Finally  when  the  estimates  7's  in  (7)  replace  7's  in  (4),  ln(3,  can  be 
predicted  in  terms  of  2,1, 

ln(j3i)  7l^»l  d"  ^m^im  (^) 


or 


=  exp(7i2n  +  ...  -t- 


(9) 
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Notice  the  difference  between  the  estimated  individual  deterioration  rate 
/?,  and  the  average  group  deterioration  rate  /3,  that  takes  into  account  tlie 
random  effects. 

In  order  to  predict  the  proportion  defective  of  a  randomly  selected  lot 
{  that  was  not  used  in  the  experiment,  one  can  use  the  following: 

p){t)  =  exp0o+/3^t)/(l  +  exp0o+0/t))  (10) 

where  =  exp(7i2,'i+, +imZem)- 

Subsequently,  one  can  estimate  the  expected  time  t)  when  the  quality 
of  ammunition  lot  H  reaches  a  predetermined  level  p: 

i\.{p)  =  [ln{p/{\  -  p))  -  3o]/3i..  (11) 

3.  STEP-STRESS  MODEL 

In  order  to  formulate  the  deterioration  model  with  step-stress,  the  fol¬ 
lowing  assumptions  are  made  based  on  Nelson  [7]:  (1)  The  remaining  life 
of  ammunition  depends  only  on  the  current  cumulative  proportion  defec¬ 
tive  and  the  current  stress  associated  with  the  depot  regardless  how  the 
proportion  is  accumulated;  (2)  If  held  at  the  current  stress,  defective  items 
will  occur  according  to  the  logistic  function  for  that  stress  but  starting  at 
the  previously  accumulated  fraction  failed;  (3)  The  change  in  stress  has  no 
effect  on  life  -  only  the  level  of  stress  does. 

Now,  suppose  that  ammunition  lot  i  would  be  stored  in  depot  ki  dur¬ 
ing  [0,  ti],  in  depot  k2  during  [<i,  ^2]  and  finally  transferred  to  depot  k^  and 
stored  during  [<2,  4]-  In  this  case,  a  routing  sequence  of  locations  of  depots 
becomes  {ki,k2,  ks)  along  with  associated  duration  [0,  4],  [4,  4],  [4,  4],  re¬ 
spectively. 
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First  of  all,  we  define  the  predicted  cumulative  proportion  defective  of 
lot  i  in  the  /th  sequence  of  storage,  depot  ki,  as  p^P(t).  For  lot  i  that  would 
have  been  stored  in  depot  k]  from  the  beginning  to  time  ti,  the  expected 
cumulative  proportion  defective  by  time  t  is  predicted  as 


.  '(fci)  .  -(ki) 

Pi{i)  =  P,-  (<)  =  exp(3o  +  4.  *)/(!  +  exp{0o  +  0,  0)  (12) 

-Aki) 

where  0  <  t  <  ti  and  /3,  =  exp(%0. 

Next,  ammunition  lot  i  is  moved  to  depot  k2  from  depot  k^.  Depot  A;., 
has  an  equivalent  starting  time  S]  which  would  have  produced  the  same 
proportion  defective  as  in  depot  ki  at  time  ti  if  the  ammunition  lot  had 
been  stored  in  depot  k^  from  the  beginning.  Such  Si  would  satisfy  the 

*  -Akj)  ,  -.(tj) 

following  relationship:  p\  \si)  =  p\  ’(t\)  or  30  +  0,  Sy  =  3o  +  3i  t:-  Thus, 

-Aki)  -Ak2) 

S.  =  (3,  /P,  )<,.  (13) 

As  a  result,  the  predicted  cumulative  proportion  defective  of  lot  i  in 
depot  Ajo  by  time  t  after  transferred  from  depot  A;i  is 


-  *(^2)  .  1(^2) 

Pi{^)  =  Pi  {t-tl+Sy)  =  exp{3o+0i  {.t-tx+Si))/{l  +  exp{3o+0,  {t-ty+S,))) 

(14) 

where  ty  <  t  < 

Similarly,  at  the  third  sequence  of  the  routing,  depot  k^  has  the  equiv¬ 
alent  starting  time  S2  which  satisfies  pf\s2)  =  P^Hh  —  ti  +  Sy).  Thus,  Sj  is 
Akz)  ,  -Ak2) 

the  solution  of /3o  + /3.  S2=0o+0i  (h  -  U  +  sy),  i.e., 

-AkA  -Aks) 

^1  =  {0,  10  i  ){h-~h+Sy).  (15) 
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The  cumulative  proportion  defective  of  lot  i  in  depot  by  time  t  after 
transferred  from  depot  k\  and  depot  k^  sequentially  is  predicted  as 

.  -Ah)  ,  Ah) 

Pi{i)  =  Pi{t-h+S2)  =  exp{f3o+(3^  {t-ti+S2))/{l+exp(J3o+0  ^  (i-ts+so))) 

(16) 

for  time  period  la  <  t< 

In  general,  for  lot  i  which  is  transferred  from  depot  ki^y  to  depot  ki  at 

time  ti_i,  depot  fc;  has  the  equivalent  starting  time  which  satisfies  the 

Ah)  .  Ah-i) 

following  relationship;  /3o  Si-i  =  /?o  (<i-i  —  ti-a  +  si-a)  where 

/  =  2, L,  A)  =  0  and  Sq  =  0.  Therefore 


Ah-i)  Ah) 

Si-1  =  (/3,  //3,  )(A-i  ~  A-2  +  S/-2)-  (17) 

Consequently  the  cumulative  proportion  defective  of  lot  i  in  depot  ki  by 

time  t,  ii-i  <  t  <  ti  is  estimated  as  follows: 

Piit)  =  Pf{i-  ti^l+Si-i) 

Ah)  *  Ah) 

=  ea:p(5o+/3.  (t- A-i +Si-i))/(l  +  ea:p(/?o +/?,  (t- A-i  +  s;-i)))  (18) 

for  time  period  A-i  <  <  <  A- 

Thus,  the  predicted  cumulative  proportion  defective  of  ammunition  lot 
i,  Pi{t),  for  step-stress  pattern  (/ci,  A:2,  fca)  along  with  eissociated  duration 
[0,  <i],  [ti,  ta],  [A>,  4]  consists  of  segments  of  the  pf  \  ),pf\  ),  and  pf\  ). 

4.  ILLUSTRATION 

In  order  to  illustrate  implementation  procedures  for  the  methods  sug¬ 
gested  in  this  paper,  a  numerical  example  is  generated  based  on  the  pa¬ 
rameters  used  in  the  guideline  of  the  U.S.  Army  Ammunition  Surveillance 
Procedure  [12], 
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Suppose  20  ammunition  lots  were  used  for  experiment  to  estimate  /3,’s 
where  (i  =  1,..,20).  The  experimental  lots  are  exposed  to  four  different 
levels  of  constant  stress  representing  the  conditions  of  depot  1  (lots  1,..,5), 
depot  2  (lots  6, ..,10),  depot  3  (lots  11,. .,15),  and  depot  4  (lots  lt),..,20).  All 
the  lots  are  inspected  annually  based  on  sample  size  of  20  (uy ).  Inspection 
starts  when  the  lot  is  3  years  old  ant  it  is  done  thereafter  every  three  years 
until  15th  year.  Table  1  contains  information  regarding  the  series  of  the 
number  of  defective  items  (y^)  found  in  20  ammunition  lots  (AT). 

This  information  is  used  to  obtain  ML  estimates  of  parameters  in  the 
within-lot  model  (1)  0o,  /3i,..,/?2o  and  A).  For  this  step,  PROC  LOGISTIC 
of  a  statistical  package  SAS  [8]  is  used.  In  Figure  1,  sample  patterns  of  the 
actual  deterioration  (ActualrXy  /  )  are  overlaid  to  those  of  the  estimated 

A  A  A  A 

deterioration  (Fitted:  (exp(/3o  +  /?i<i>)/(l  +  exp(/3o  +/3,ty)))  against  time 

t,)- 

The  between-lot  model  (2)  is  formed  to  relate  deterioration  rates  to  the 
depot  characteristics,  ln/3,  is  used  as  dependent  variable  and  the  four  dum¬ 
my  variables  (2,i,..,2ri)  representing  4  different  depots  are  used  as  covariates 
without  an  intercept:  z,i  =  1  for  ammunition  lots  stored  in  depot  1,  oth¬ 
erwise  Zji  =  0;  Zi2  =  1  for  ammunition  lots  stored  in  depot  2,  otherwise 
Zi2  =  0;  Zq  =  1  for  ammunition  lots  stored  in  depot  3,  otherwise  Zq  =  0; 
and  Zn  =  1  for  ammunition  lots  stored  in  depot  4,  otherwise  z^  =  0.  Given 
such  (z,i,  ..,zn)  as  well  as  (/?',  and  A),  the  IMSL  subroutine  ZSPOW  [4],  is 
applied  to  (6,7)  in  order  to  obtain  71, ..,  74  and  {Z'  WZ)~^.  These  estimates 
are  summarized  in  Table  2. 

Now  we  apply  these  results  to  predict  the  proportion  defective  of  an 
ammunition  lot  i  which  would  be  transferred  from  depot  1  to  depot  4 
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following  the  route  given  in  Table  3.  The  route  given  in  Table  3  can  be 
related  to  the  following:  ammunition  lot  i  would  be  stored  in  the  temporary 
depot  1  for  1  year  and  transferred  to  the  permanent  depot  2  where  it  would 
be  kept  for  the  next  4  years.  The  lot  i  then  would  be  sent  to  the  depot  3  in 
warship  and  stored  for  1  year  and  brought  back  to  depot  4  where  it  would 
remain  until  ultimate  usage. 

The  expected  proportion  defective  in  depot  1  by  time  t,  0  <  t  <  1,  can 
be  predicted  as 


P.(<)  =  P^\^)  = 


ea:p(— 5.8551  +  e2:p(— 1.6058)  t) 


(19) 


(1  +  ea:p(-5.8551  +  ea:p(-1.6058)t))  ‘ 

The  expected  proportion  defective  of  ammunition  lot  i  stored  in  depot 
2  during  1  <  <  <  5  after  being  transferred  from  depot  1  is  predicted  as 

f,H\  +  exp(-1.5432)(t-  0.0606)) 

^  ^  ^  (l  +  eip(-5.8551  +  eip(-1.5432)((-  0.0606))) 

(20) 

Next,  the  expected  proportion  defective  of  ammunition  lot  i  stored  in 
depot  3  during  5  <  i  <  6  after  depot  1  and  depot  2  is  predicted  as 

t,  -  n  oniB^  -  +  exK-1.3566)(t-  0.9016)) 

^  ^  (H-exp(-5.8551  +  exp(-1.3566)(t- 0.9016)))' 

(21) 

Similarly,  the  expected  proportion  defective  of  ammunition  lot  i  stored 
in  depot  4  during  t  >  6  after  storage  in  depot  1,  2  and  3  is  predicted  as 


2.5887)  = 


exp(-5.8551  +  ea:p(-0.9547)(t-  2.5887)) 

(1  +  exp(-5.8551  +  ea:p(-0.9547)(t  -  2.5887))) ' 

(22) 


Figure  2  shows  the  segments  of  the  estimated  cumulative  proportion 


defectives  obtained  under  the  four  different  depot  conditions.  First,  the 
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quality  of  the  ammunition  lot  deteriorates  according  to  the  average  condi¬ 
tion  of  depot  1  up  to  time  t\.  After  it  is  transferred  to  depot  2,  the  stockpile 
deteriorates  following  the  pattern  fitted  for  depot  2,  starting  at  the  accu¬ 
mulated  proportion  defective  due  to  the  condition  of  depot  1.  Similarly, 
when  the  lot  is  shipped  fi'om  depot  2  to  depot  3,  from  depot  3  to  depot  4, 
deteriorating  patterns  follow  the  corresponding  logistic  models  fitted  under 
each  depot  condition,  respectively  starting  at  the  accumulated  proportion 
defective  at  the  previous  depots.  Figure  3  gives  their  connection  represent¬ 
ing  the  predicted  cumulative  proportion  defective  of  the  ammunition  lot 
under  the  step-stress  described  in  Table  3. 

Under  the  route  given  in  Table  3,  the  proportion  defectives  of  ammuni¬ 
tion  lot  i  in  depot  4  by  t  =  10  is  predicted  as  follows: 

P.(IO)  =  #(7.4113)  =  0.0473.  (23) 

Similar  calculation  provides  PiO.2)  =  0.0969,  Pt{lS)  =  0.1365,  and 
p.(14)  =  0.1880. 

When  the  substandard  quality  of  the  ammunition  lot  is  set  at  p  =  0.15, 
this  ammunition  lot  should  be  used  or  replaced,  at  the  latest,  by  the  end 
of  13th  year  of  storage  in  depot  4  provided  it  had  followed  the  route  given 
in  Table  3. 
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5.  DISCUSSION 

A  two-stage  random  effect  logistic  regression  analysis  is  applied  to  pre¬ 
dict  the  declining  quality  of  ammunition  stockpile  when  a  series  of  location 
and  the  duration  of  storage  associated  with  the  ammunition  lot  is  given. 
First  of  all,  average  deterioration  rates  are  estimated  using  constant  stress 
model.  Next  the  prediction  method  for  the  cumulative  proportion  defective 
of  an  ammunition  lot  is  described  using  a  step-stress  model. 

We  use  a  logistic  regression  model  to  fit  deteriorating  patterns  of  ammu¬ 
nition  stockpile.  The  patterns  of  dechning  qualities  are  observed  based  on 
sampling  inspection.  Other  alternatives  to  the  logistic  model  include  com¬ 
plementary  log- log  model  and  probit  model  (McCullagh  [6]).  For  instance, 
the  complementary  log-log  and  the  probit  within-lot  models  which  corre¬ 
spond  to  the  logistic  within-lot  model  (1)  would  be  Py  =  1  —  ea:p{  —  exp{(3o  + 
0ity)),  and  Pij  =  ^(/3o  + respectively,  where  ^(.)  is  the  standard  nor¬ 
mal  integral.  Even  in  these  model  specifications,  parameters  /3,’s  represent 
the  deterioration  rates  as  described  in  the  logistic  regression  model  (1).  To 
estimate  /?,’s,  the  maximum  likelihood  estimation  methods  given  in  (3)  can 
be  applied.  The  corresponding  between-lot  model  analyses  are  essentially 
the  same  as  described  in  the  two  stage  estimation.  Choice  of  the  specific 
within-lot  model  among  logistic,  probit  and  complementary  log-log  models 
depends  on  the  pattern  of  data  gathered  fi’om  the  experiment. 

When  the  response  is  taken  in  terms  of  the  change  in  a  certain  attribute 
of  an  item  such  as  water  content  in  propellant,  rather  than  counting  the 
number  of  defectives  based  on  a  certain  sampling  scheme,  a  nonlinear  mod¬ 
el  can  be  used  to  formulate  deteriorating  patterns  of  ammunition  lots.  A 
typical  example  is  use  of  a  negative  exponential  growth  curve  model.  How- 


14 


ever,  the  selection  of  appropriate  nonlinear  model  again  depends  on  the 
observed  declining  characteristic  of  the  quality  of  ammunition  stockpile. 

In  this  paper,  different  levels  of  depot  conditions  are  viewed  step-stress. 
Degradation  is  assumed  to  be  independent  of  the  sequence  of  step-stress, 
(i.e.,  sequence  of  depot  routing).  Sometimes  abrupt  changes  between  two 
consecutive  depots  may  cause  faster  deterioration  of  stockpile  than  expect¬ 
ed  in  these  two  depots.  Developing  models  that  accommodate  sequential 
dependence  is  left  as  one  of  areas  for  further  study. 
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Table  1;  Number  of  Defectives  i/y  Observed  in  Sample  Size  of  20  at  the  jth 
Inspection  of  Lot  i _ 


group  1 

lot  i 

group  2 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

U 

12 

0 

1 

1 

0 

0 

1 

0 

1 

1 

1 

15 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0, 

0.1494 

0.1966 

0.1966 

0.1494 

0.1494 

0.1966 

0.1494 

0.1966 

0.1966 

0.1966 

se{/30 

0.0813 

0.0644 

0.0644 

0.0813 

0.0813 

0.0644 

0.0813 

0.0644 

0.0644 

0.0644 

group  3 

lot  i 

group  4 

iv 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

1 

0 

1 

1 

1 

1 

9 

1 

0 

0 

0 

2 

1 

2 

2 

3 

2 

12 

1 

1 

1 

1 

3 

3 

4 

4 

5 

4 

15 

1 

1 

1 

1 

5 

5 

6 

6 

8 

6 

0i 

0.2197 

0.1966 

0.1966 

0.1966 

0.2836 

0.3827 

0.3520 

0.3520 

0.3&)27 

0.3520 

se0i) 

0.0583 

0.0644 

0.0644 

0.0644 

0.0471 

0.0414 

0.0421 

0.0421 

0.0414 

0.0421 
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Table  2:  Fitted  Between-Lot  Model 


1 

Z'fVZ 

1 

-1.6058 

0.0631 

0.0316 

0.0257 

0.0163 

2 

-1.5432 

0.0316 

0.0490 

0.0236 

0.0149 

3 

-1.3566 

0.0257 

0.0236 

0.0297 

0.0121 

4 

-0.9547 

0.0163 

0.0149 

0.0121 

0.0106 

Table  3:  Scenario  for  Location  and  Duration  of  Storage 
Depot  12  3  4 

Period  [0,1]  [1,5]  [5,6]  [6,—] 
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Cumulative  Proportion  Defective 


Figure  2.  Segments  of  the  Predicted  Cumulative  Proportion 
Defective  of  the  Ammunition  Lot  under  Constant  Stress 
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